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ABSTRACT 


A simple means for bounding the propagated trajectory 
error induced by an impulsive initial error in the n~body 
field is developed. Formulas are derived for the maximum 
acceptable initial error, given the mission tolerances, and 
for the maximum time of negligible propagated error, given 
the initial impulses and mission tolerances. An example 
illustrates the use of the results within the solar system. 
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EXPLICIT UPPER BOUNDS FOR ERROR PROPAGATION 
IN THE n-BODY FIELD 


by 

Albert Rosen* 

Goddard Space Flight Center 


INTRODUCTION 

One of the central problems of astronautical guidance is to estimate what effect a change in 
the position and velocity of a spacecraft at some initial time has on its position and velocity at 
some later time. Such a change may occur accidentally, as when there is an error in the value 
assigned to the initial state vector, or on purpose, as when a rocket motor induces an impulsive 
velocity change. 

If the spacecraft has an essentially Kepler orbit, the propagation of the initial error can be 
analyzed precisely by matrix methods (References 1 and 2). If this error arises from a rocket 
motor, taking the cumulative velocity change as an impulse will cause the solution to be approx- 
imate. The closeness of approximation is discussed by Robbins (Reference 3). In circumstances 
where the Keplerian approximation is too crude, it is still possible to use matrix methods for the 
expression of error propagation, provided that the error is small enough to permit valid lineari- 
zation of the relative force between the deviated and initial trajectory. Battin (Reference 4) and 
Baker (Reference 5) develop the proper forms for the matrices involved, and Dusek (Reference 6) 
gives a heuristic discussion of the validity of such linearization procedures. 

Let the propagated error vectors be defined as the differences of the position and velocity 
from their values on the original orbit, at some time after the initial error is applied. Then, if 
it is desired to represent these vectors accurately, one of the above methods may be used. 

It may be that certain errors are so small that they will produce no appreciable change on the 
trajectory, even after a considerable time. One way to verify this would be to include the effects 
in a numerical integration scheme and compare the resultant trajectory with a trajectory com- 
puted without considering the effects. Another way would be to use the linearized equations of 
motion in a matrix propagation. But it would be far simpler if there were some easy way to place 
an upper bound on the propagated effect of such a small error; then it should be enough to note that 
this upper bound is small (in view of the required precision). 


♦Student: St. Josephs College, Pennsylvania. 
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We can develop a simple expression for such an upper bound by assuming that 


1. The neglected error is approximated by the instantaneous variations r 0 and v Q in position 
and velocity at time t = 0. 

2. The motion of the spacecraft is adequately described by the n-body equations, for the time 
of interest: 

3. Enough is known about the initial or "unperturbed" trajectory to permit placing lower 
bounds on the closest approach distances of this trajectory to each of the n bodies, in the 
time interval under consideration. 

As a measure of the effect of the error induced by the initial error r 0 , v 0 , in position and ve- 
locity, we define two vectors as follows: r(t) is the difference between positions at time t, in the 
perturbed and unperturbed orbits, respectively. v(t) is the same difference for velocities. At 
each time t, the vector Q. will be the instantaneous vector to the i th body in the system, from the 
unperturbed trajectory. 


AN INEQUALITY FROM THE EQUATIONS OF MOTION 

Using the above definitions, Newton's Law of Gravitation leads at once to: 



7(0) = 7 0 ■ 

v(0) = v 0 , (1) 


where 

fx i = gravitational constant of 1 th body, 
n = number of bodies in the system. 


Eventually we will bound the magnitude of the solution of this highly nonlinear equation by the solu- 
tion to a certain linear homogeneous scalar differential equation with constant coefficients, which 
can be derived from it. 

Assume that | r(0) | < |q. (0)1 (i = 1, 2, • • • , n). Then, for a short enough time, I r(t)| < |q. (t)| , 
because all these functions are continuous in time. 
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We first show that | r(t)| < | (t)| implies the inequality 


Q ; +? Qi 

1 

lQi+7| 3 iQi 1 3 

1 (|QJ -I?l) 


( 2 ) 


To verify this, expand the square of the quantity on the left, denoted by U, in powers of |Qj \/\p\ and 
Irj/lQ. |, where ~p = Q. + r. The result is 


U 2 



i iQjl 1 iQjl 4 l iQil 3 
4 lp| + 4 l?l 4 " 4 l?l 3 



The maximum value of U 2 occurs when 


S( x ) 




a 


(3) 


maximizes, x - \Q i \/\p\ 7 a □ l-|r|^/|Q. I 2 and the magnitudes |Qj and I r | are fixed. Since 
I r | < |q. | , the inequality 4 restricts a and x: 


0 < a < i and 


1 

— - < X < 

i + VTi -«) - - 


i 

i-yrr^o ’ 


(4) 


with the extreme values of x in inequality 4 corresponding to r parallel and antiparallel to Qj . 
But the derivative of S(x) is 


S' (x) 



(5) 


If we equate S' (x) to zero and solve for x, the discriminant is (a 3 +4)/256, which must be posi- 
tive, according to inequality 4. Hence S' (x) = 0 has one real root. Evaluating S' (x) for the 
extreme values of x in inequality 4 gives 


and 



il a (l + 2) 


which is negative, 



which is positive 


( 6 ) 


so that the single real root for S' (x) =0 lies inside the interval of x given by inequality 4. 
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This means that S(x) takes on local maximum values for x - l/(l+yi - a ) and x = l/(l -/l - a ) 
because, where S' (x) = 0, we have an absolute minimum for S(x). Consider that S(x)-*oo as x-»+oo. 

Hence, either for r parallel to Q A (giving U = U p ) or for r antiparallel to Q. (giving u = U a ), U 
takes on its maximum value. Evaluating u for these two orientations of r and Q. gives 


(IqJ - I r|) 2 I Q ± 1 2 IqJ 2 |_ IqJ IqJ 


r r 

2 -zr 1 - + 3 


... (n + 1) 


Qi I n 


(7) 


(l Qi I ~\r\) 2 IqJ 2 IqJ 2 L IqJ Iq, 


' 3 


(-l) (n+1) (n+ 1) 4r~ + **• * (8) 

I Qi I n 


so that U a >U p , and U a is an upper bound for U for all orientations of T and Q. . This completes 
the proof of inequality 2. 

Returning to Equation 1, the triangle inequality can be combined with inequality 2 to give the 
result 


n 

± L 


V-x 


(|Q,[-IT[) : 



provided that |r| < IqJ . This is identical with the inequality 

|7| £ |7| f 2 (l?l. IQJ) , 


where f is a function of the n + 1 variables x, , y 2 , * • • , y n given by 


f(x, y t ) 



(9) 


(10) 


(ID 


It is assumed in Equation 11 that 0 <x<y £ (i = 1, 2, • • • , n). Clearly f and f 2 are monotonically 
increasing functions of x. Also, f df/dx is a positive monotonically increasing function of x, be- 
cause it is given by 


d f 


f ( x ’ y 0 ^ ( x ’ y i) = 


^ (>-£)' 


( 12 ) 
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Both f (x, y.) and f (x, y £ ) <9f (x, y.^x are monotonically decreasing in all the y^s, so long as x 
remains positive and less than each of these variables. 

For sufficiently small t, inequality 10 is valid and can be combined with 

|r| < [ f |r(s)| ds 2 + |rj + | v Q | ■ t , (13) 

Jo Jo 


which is true for all t , to give: 


|r| < f f I r(s) | f* (|r|, |Qj) 
Jo Jo 


ds 2 


(14) 


This implicit integral equality enables us to get an explicit upper bound on |r| and I v| . 

AN EXPLICIT UPPER BOUND FOR |r {AND | v| 

Consider the differential equations 

w - 
w(0) = 

w(0) = 

For the present, x is a parameter which has the range 0 < x < . n , where C mi n is the smallest of the 

C ,'s. We assume that the C/s are chosen on some interval 0 < t <Tj in such a way that 

C; < | Q ; (t)| ; 0 < t < T, , (16) 


V f2(x, Cl ) 


(15) 


and we restrict r 0 to satisfy the relation 


0 < 


< c; 


(17) 


This last restriction causes no practical difficulties, because the C { represent lower bounds on the 
closest approach of the unperturbed trajectory to the i th body, and no initial error in position will 
in practice be as large as any of these distances. 

The solution of Equation 15 on 0£t <Tj is given by 

|v Q | sinh f (x, C.) t 

w = | r Q | cosh f (x, Cj t + f ( x , C.) • ( 18 ) 
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which satisfies the integral equation 


w( x, t) 


n w(s) f 2 (x, Cj ds 2 + |rj + |v 0 | t 

' 


(19) 


0 < t < T 


Comparing Equation 19 and inequality 14 leads us to suppose that w(x, t) may provide an upper 
bound for |r(t)| , provided that x can be chosen as an upper bound on I r(t)| , during the time under 
consideration. If this were so, we might try to choose x to satisfy 


w(x, T) < x ; 


T < T, , 


( 20 ) 


which would be consistent, at least, with the supposition that w(x, t) provides an upper bound on 
| r| on 0 < t <T, T < Tj, and that x provides an upper bound on I ?l on that interval. 

These heuristic considerations lead to the following theorem. 

Theorem I. Let w(x, t) be defined as in Equation 18, let the C/s be chosen to satisfy Equation 16, 
and let | 7J be restricted as in inequality 17. Then if x satisfies inequality 20: w(x, t) provides 
an upper bound on |r| on 0< t <T, and 3w(x, t)/dt provides an upper bound on |v| on that interval. 

We show that w (x, t) provides an upper bound on I r| on the half-open interval 0 < t < T. The 
continuity of |r| and w(x, t) in time then enables the extension of this result to the endpoint of that 
interval. 

Assume the contrary, that there is a time in 0 < t <T such that w(x, t) < | r(t)| . Then, by the 
continuity of | r| and w in time, there is an open subinterval of Of t <T where w(x, t) < 7 ( t ) | , and 
there is a first such subinterval. Letting denote the beginning of the subinterval, y(t) = | r(t)| 


time 

> tj, sufficiently close to t p 

(a) 

7(t) < 0 

0 

1 A 

rf 

|A 

r+ 

(b) 

y(t) = 0 

rh 

II 

rf 

(c) 

r(t) > o 

‘,<‘1*2 

(d) 

Hi 

! A 
X 

0< t < tj 

(e) 

r(t 2 ) > y(t) 

t, < t < tj 


(f) 


t, < t 2 < T. 


(21) 


To see that inequality 21d is true, note that inequality 21b implies that | r (tj| < x, because 
w(x, t,) <w(x, T)<x. Continuity then demands that I r(t)| remain less than or equal to x for a 
short enough time after t 2 — i.e., for a t 2 close enough to t,. 
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I 


Inequality 21e follows from the fact that 


r(t) 


I 7(s)ds, 


tj £ t < t 


2 • 


since y (t J = 0. But since this is positive for every tin t 1 <t<t 2? then y( s) must be positive in 
that interval, so y(t) is monotonic ally increasing in t in that interval. 

Since inequality 21d is true, and x£ |q. (t)| onO<t<T,we conclude that inequality 14 must 
hold in 0 £ t < t 2 . Subtracting Equation 19 from inequality 14 then gives, letting t = t 2 , 


0 < 


y(t 2 ) < | | [l r(s)| f 2 (l r(s)| , |Qi (S)|) - w(x, s) f 2 (x, C.) 


ds 2 


( 22 ) 


since the contribution to the integral from 0< t < 1 1 is not positive. But by inequality 21d, and the 
fact that C i < |Q. | (i = 1,2,..*, n), we can reduce this to 

o < y(t 2 ) £ f 2 (x, Cj) [ f [|r(s)| -w(x, s)j ds 2 . (23) 

Jt i 

The mean- value theorem then gives, with 1 1 £ t 3 £ 1 2 ; 

o < r(t 2 ) < f 2 (x, c.) • r(t 3 ) (t 2 ~ t i) 2 • (24) 

But by inequality 21 e, this implies that 

o < r(t 2 ) < f 2 (x, Cj) y(t 2 ) (t 2 - t,) 2 . (25) 

This is a contradiction, as may be seen by dividing through by y(t 2 ) and letting (t 2 - approach 
zero. Therefore w(t ) > | r| on 0£ t < T and, by continuity, onO£t£T. 

To see that dvj/dt > |v| on 0 < t < T, we use this result, noting that f 2 (|r| , |q. |] £ f 2 (x, c.) on 
0 £ t £ T, and also that inequality 16 is true on 0 £ t < T. But, always: 



(26) 


II 
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Taking the first partial derivative of Equation 19 with respect to time, and combining this with 
inequalities 26 and 10 gives 


I v[ - ^-^Tt t) < J f2 (l r| , I Qi l) " Wf 2 (x, C.)J ds 
which by the above remarks is reducible to 

I v| - ^ £ f2 ( x > c i) [l r(s)| - w(x, s)J ds . 


( 27 ) 


(28) 


Since the integrand has been shown to be non- positive for all t in 0< t < T, the left-hand side of 
inequality 28 must also be non-positive in that interval. This concludes the proof of theorem I. 

The question now is: for what range of values of T, | r 0 | , and | v 0 | does there exist a range 
of x values that satisfy inequality 20, and what is this range of x values? Can these x-values be 
conveniently found? 

Consider the following expressions for w(x, T) and <9w(x, T)/dx: 


■ _» t _ sinh fT 

w(x, T) - | r Q | cosh fT + v Q f 


= Ir 0 l + |v 0 | T + 0 (T 2 ) 


<9w(x, T) 
dx 


, 0 /sinh fT\ - - _ If sinh fT\ 

I r 0 l T 2 (ff ) fT j + I v 0 | T 3 (ff ) ^^7 ^cosh fT - —Jj J 


{• 


(ft ) 2 (fry 4 


|r 0 | T 2 ( f f ' ) <1 + ^ff~ + ^- 5 f- + 


(29) 


+ |v 0 | t 3 ( f f ' ) [(^T - jr) + ( fT > 2 (tT - ^f) + (fT ) 4 (jr- yr) 


(30) 


In Equations 29 and 30, f - f (x, and f' - <5f(x, C-^/dx. 

Both series in Equation 30 are monotonically increasing in fT, and both f and ff' are mono- 
tonically increasing in x; therefore <9w(x, T)/dx is monotonically increasing in x and T. 

When x = | r 0 |, w(x, T) - x is positive for T > 0, and w (x, T) - x is also positive for x-c min , be- 
cause this makes f (x, C.) infinite. The fact that < 3 w(x, T)/dx is positive and monotonically increas- 
ing in x implies that the function w(x, T) - x is always curving upward in x, so that it equals zero at 
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none, one, or two points. If there are any roots, then w(x, T) - x has a minimum which is also an 
extremum on | r 0 | j< x < C min . Clearly then, it is possible to find x such that w(x, T) < x is satisfied 
if and only if there are real roots to 


w(x, T) - X = 0 ; |r 0 | < x < C min , 


(31) 


and any value of x between or including these roots will satisfy w(x, T) < x. 

For small T, Equation 31 has two roots, since the series in Equation 29 can be made arbitrarily 
close to | r 0 1 by choosing a small enough T, regardless of where x is in I r 0 1 < x < C mi n . For a large 
T, the value of <9w(x, C^/dx at x = | r 0 | becomes large, and since it is always increasing in x on 
I r J < x < c . , there can be no intersection points of w(x, T) with x; and no roots to w(x, T) - x = 0. 
There is a critical value of T, therefore, which permits one and only one solution to Equation 31, 
and this is the largest T for which inequality 20 can be satisfied. At this critical value of T, the 
root of Equation 31 defines the extremum of w(x, T) - x 3 so the condition 


JZ [ w (x, T) - x] = 0 , 


x < C 


(32) 


must also be satisfied. Solving Equations 31 and 32 simultaneously for T gives the desired critical 
value, provided that we know I v 0 | , |r 0 | , and all the C i T s. 

A simpler method for finding the critical T can be derived, however. The function w(x, T) has 
the form w = w(fT, f). The total derivative of Equation 31 with respect to x is therefore given by 


dw(x, T) dT 1 <9w( x, T) . dj_ 

d( fT) dx ^ dx * d f (9x 


(33) 


where we use, the fact that Equation 31 determines T as a function of x. On the other hand, Equa- 
tion 32 is identical with 


<9w r<9 f 1 dw d f 

d( fT) L<9 x ‘ T J + df dx " 


(34) 


Since there is a unique solution to Equations 31 and 32 on |r 0 | < x <C nin , and since f • (<9w/d(fT)) 
is readily verified to be nonzero, we conclude that Equation 31 combined with Equation 32 is the 
same as Equation 31 combined with the condition that dT/dx = 0. This condition determines the 
maximum value of T(x) on |r 0 l £ x <c min > because T must vanish at both endpoints, in order that 
Equation 31 shall have roots. Solving Equation 31 for T explicitly leads to the following theorem 
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Theorem 77. Let t(x) be defined by 


T ( x ) = f (x, Cj) ^ 


r 

1 V 0 1 1 

2 - i r oi + 

f (x, C.)_ 


> • 


(35) 


f(x, C.) 




Then, for all T satisfying the inequality 

T < max(r(x)') , |r Q l £ x < C min , (36) 

there are values of x satisfying w(x, T) < x; all such values lie in the interval x 1 < x < x 2 , where x x 
and x 2 are the two roots of the transcendental equation: 

|v 0 | , 

| r 0 ! cosh f (x, C^) T + f c ^ sinh f (x, C.) T = x , (37) 


on I r o I £ X < C m ; n . 

Suppose that the impulse v 0 and position error r 0 were assigned and it was required to know 
for how. long a time the propagated effect of these vectors would remain less than x in position mag- 
nitude (| r 0 | £x<C min ) . The answer to this is: set T 11 r(x). This is legitimate, because, if 
T = r(x) , then Equation 37 is satisfied, so w(x, T) = x. Hence inequality 20 is satisfied, and W(x, t) 
bounds the error in position on 0 < t < T, and x bounds the error in position on 0 £ t < T also, as 
needed, since x>w(x, t) onOft <T. 

Another theorem may be presented here: 

Theorem III . If it is required that the propagated position error be less than or equal to x, where 
|r | < x <c min , then this is assured by taking the time of propagation equal to r(x), where r(x) is 
defined in Equation 35. 

Suppose, next, that it is desired to know how large an initial error, r 0 and v 0 is permissible, 
in order that the propagated effects of these errors over a pre- assigned time T should be less than 
*(l?„l <*<C mln ). 

To answer this, note that for any x satisfying 0£x <C min , it is possible to choose the coef- 
ficients of cosh f (x, T and J^sinh f(x, C £ ) (x, cj so that 

. V / sinh f (x, C.) t\ 

|rj (coshf(x, Cj) • Tj + |v 0 | [ " f j < x (38) 

is satisfied, thereby guaranteeing that x is indeed an upper bound on I r| after time T. In particular, 
if | r 0 1 = 0, then the propagated error is less than x after time T, provided that I v 0 1 satisfies the 
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inequality 


_ , xf ( x - c i) _ 

V °' - sinh (f (x, C.) t) 


( 39 ) 


Note that although r(x) in Equation 35 may become formally greater than T u the time for which 
the C.fs were determined, the fact has no physical significance. If this happens, we can only be as- 
sured that x exists to satisfy inequality 20 on 0 £ t < Tj . If it is desired to extend the interval past 
Tj, then it is necessary to find Ci 's for this extended interval to make all the results valid. 

It is also worth noting that looser tolerances in mission requirements (larger x) in Theorem 
IE do not necessarily produce larger times of negligible effect, because t(x) does have a maxi- 
mum. If x is chosen to be larger than the value which maximizes ^(x), then a looser tolerance will 
decrease r(x), giving pessimistic results. In such a case, the solution is to take the time as max- 
imum r(x), since the error x corresponding to this maximum is certainly less than an error cor- 
responding to the times less than the maximum, brought on by large x. 

Finally, if x/C min is small, then, to a good approximation, f (x, = f(o, C^) , i.e., 



Example 


A spacecraft enroute to the outer planets passes through the asteroid belt and may be subjected 
to an impulsive velocity change there as large as 50 km/day = 5.8 x 10~ 4 km/sec. This figure may 
have been arrived at by placing lower bounds on the closest approach of the spacecraft to the as- 
teroids (as determined by their physical radius), lower bounds on the relative velocity of rocket 
and asteroid, and upper bounds on the masses of asteroids, and by assuming Keplerian interactions 
for the short encounter time with the asteroid. We wish to examine the orbit for some part of its 
path during which the closest approach distances to the sun and other large bodies in the solar sys- 
tem is given by Table 1. 


The other bodies in the solar system may be considered negligible in influencing the propaga- 
tion of the impulse. For how long does the propagated error in position remain less than 0.008 AU? 


Solution: Use Theorem HI. First calculate f (x, c^. Since 0.008/C i « 1 by Table 1 for all the 
bodies, we use the approximation: 


f (x, Ci 



where t± s = Sun T s gravitational parameter. 



(41) 


11 



Table 1 


Closest Approach of Spacecraft to 
Bodies in Solar System 


Body 

Closest Approach, AU* 

Sun 

1.50 

Mercury 

1.00 

Venus 

1.60 

Earth 

1.00 

Mars 

0.80 

Jupiter 

3.50 

Saturn 

8.00 

Uranus 

9.00 

Neptune 

14.00 


*1 AU = 1.49 x 10 8 km. 


Now the largest value of Mi/^ s is that for 
Jupiter, which is 0.9 x 10“ 3 . Also, by Table 1, 
the largest value of (c s /c i ) 3 is that for Mars, 
which is 6.6. Hence, to a percent accuracy, 


f(x, Ci ) 



where m s is in AU 3 days 2 , C s is in AU. That 
is: 


f(x, C.) 


/ 


5.92 x 10~ 4 
( 1 . 50 ) 3 


- 4475 days ' 1 . 

Therefore we are sure of a propagated error less than 0.008 (or 1.20 x 10 6 km) for a length of time 
given by 


T = r (l . 20 x 10 6 ) 


- 44 . 5 -fn 


i. 20x10® + y(l.20xl0 ®) 2 + (50.0 x 44.5) : 
50.0x44.5 


q /l.20xl0®\ 

- 44. 5 'tn l25.0 x 44.5j 

= 44. 5 'tn ( 1080) = 44.5 x 6.98 = 311 days . 


As another example, suppose we wish to know for how large a velocity impulse the propagated ef- 
fect in position would be less than 0.008 AU over a period of 100 days, assuming that Table 1 is 
valid. 

Solution: A gain , f( x , C.) = 1/44.5 and we can use inequality 39 to obtain the upper bound on | v 0 |: 

. 1. 2 x 10® 

V ° - / 100 \ 

44. 5 sinh I 44 5 j 

_ 1, 2 x 10® 

" 44.5 sinh (2. 24) 


1. 2 x 10® 

” 44.5 • 4.64 
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1.20 xlO 6 
206 


= 5811 ^ 

= 067 km/sec. 

So impulses as large as 0.06 km/sec could be tolerated without violating mission requirements. 
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Appendix A 


Symbol List 


r 0 , v 0 : Initial errors in position and velocity. 

r(t ), v( t ) : Instantaneous errors in position and velocity. 

Q. (t): Vector from unperturbed orbit to i th body. 

. Mi : Gravitational parameter of i th body. 



■ Mr . ■ 

C. : Upper bound on | Q A ( t ) | on 0 £ t < T t . 


• , n 


C . : Smallest of the C.'s. 

min 1 


w(x, t) = |r 0 | cosh f (x, C.) t + sinh f (x, C { ) t . 

r(t) = |r(t)l - w(x, t) . 


x - 1 /x 


; -|7 0 l 


+ (U x - C i)) : 


r ( x ) = f(x, Cj) 


'in 


-y • 


I r o I + f (x, C 4 ) 


1 AU = 149 X 10 6 km. 


NASA- Langley, 1968 30 
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